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Abstract 



Employing profits data of Japanese firms in 2003-2005, we kincmatically exliibit tlie static 
log-normal distribution in the middle scale region. In the derivation, a Non-Gibrat's law un- 
der the detailed balance is adopted together with following two approximations. Firstly, the 
probability density function of profits growth rate is described as a tent-shaped exponential 
function. Secondly, the value of the origin of the growth rate distribution divided into bins 
is constant. The derivation is confirmed in the database consistently. 

This static procedure is applied to a quasi-static system. We dynamically describe a 
quasi-static log-normal distribution in the middle scale region. In the derivation, a Non- 
Gibrat's law under the detailed quasi-balance is adopted together with two approximations 
confirmed in the static system. The resultant distribution is power-law with varying Pareto 
index in the large scale region and the quasi-static log-normal distribution in the middle 
scale region. In the distribution, not only the change of Pareto index but also the change of 
the variance of the log-normal distribution depends on the parameter of the detailed quasi- 
balance. As a result, Pareto index and the variance of the log-normal distribution are related 
to each other. 

PACS code : 89.65. Gh 

1 Introduction 

Log-normal distributions are frequently observed not only in natural phenomena but also in 
social phenomena. For representative example, the probability density function P{x) of personal 
income or firm size x is considered to obey the log-normal distribution [l]-[2] 



in the middle scale region. ^ Here x is a mean value and o"^ is a variance. 

*e-mail address: ishikawa@kanazawa-gu.ac.jp 

^In the large scale region, personal income distributions and firm size those follow power-law [3]. We discuss 
this in Section 2. In the low scale region, several distributions are proposed (see Refs. [4]-[5] for instance). We 
do not discuss them in this study. 





The simplest model which describes the log-normal distribution is the pure multiplicative 
stochastic process defined by 

x{t + 1) = Rit) x{t) , (2) 

where R{t) is a positive random variable. By applying this process iteratively, we obtain 

x{t) = R{t -1) R{t-2) ■■■ R{0) x{0) . (3) 

The logarithm of this equation is 

log x{t) = log R{t - 1) + log - 2) + • • • + log R{0) + log a;(0) . (4) 

If loga;(0) is negligible compare to loga;(i) in the limit t ^ oo and 

log R{i) {i = 0,1, - ■ ■ ,t—l) are independent probability variables, logx(t) follows the normal 
distribution in the limit. As a result, this model derives the stationary log- normal distribution 
PLN(a;)- 

Eq. (2) means that the distribution of the growth rate does not depend on x{t). This is 
known as Gibrat's law [1] that the conditional probability density function of the growth rate 
Q{R\xi) is independent of the initial value xi: 

Q{R\xi) = Q{R) . (5) 

Here xi and X2 are two successive incomes, assets, sales, profits, the numbers of employees and 
so forth. The growth rate R is defined as R = X2/X1 and Q{R\xi) as 

n(T>\ \ Pir{xi,R) .f.. 
^^^1"^) = P{x,) 

by using the probability density function P{xi) and the joint probability density function 

PiRixi,R). 

As far as firm sizes in the middle scale region, however, it is reported that the growth 
rate distributions do not follow the Gibrat's law (5) (sec Refs. [6]-[8] for instance). The pure 
multiplicative stochastic process model (2) cannot be applied for explaining the log-normal 
distribution in the middle scale region. Instead, we have shown that the log-normal distribution 
can be derived by using no model such as the pure multiplicative stochastic process [9]. In the 
derivation, two laws are employed which are observed in profits data of Japanese firms. One 
is the law of detailed balance which is observed in a stable economy [10]. The other is a Non- 
Gibrat's law which describes a statistical dependence in the growth rate of the past value [9]. 

In Ref. [9], we kinematically derived the static log-normal distribution under the detailed 
balance, and analyzed empirical data in a single term 2003-2004 to confirm the derivation. In 
this study, we show data analyses not only in 2003-2004 but also in 2004-2005 and 2003-2005 to 
confirm the derivation more firmly. After that, we dynamically derive the log-normal distribution 
by replacing the detailed balance with the detailed quasi-balance proposed in Ref. [11]. By this 
procedure, the log-normal distribution is described as quasi-static. 
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2 Static log-normal distribution under the detailed balance 

In this section, we briefly review the study in Ref. [9]. After that, we empirically conflrm 
the analytic result by added data analyses. We employ proflts data of Japanese firms in 2003, 
2004 and 2005 which are available on the database "CD Eyes 50" pubhshed in 2005 and 2006 
by TOKYO SHOKO RESEARCH, LTD. [12]. 

Figure 1 shows the joint probability density function Pi2{xi,X2) of all firms in the database, 
the profits of which in 2003 (a;i)-2004 (^2) exceeded 0, xi > and X2 > 0. The number of 
the firms is "227,132". Similarly, figures 2-3 show the joint probability density functions of 
the profits in 2004 (a;i)-2005 {X2) and 2003 (xi)-2005 (^2), and the number of the firms is 
"232,497" and "197,867" , respectively. Prom Figs. 1-3, we approximately confirm the detailed 
balance which is time-reversal symmetry (xi ^ X2) of Pi2{xi,X2) [10]: 

Pl2{xi,X2) = Pl2{x2,Xl) . (7) 

Figure 4 shows probability density functions of profits in 2003-2005. The distributions are 
almost stable and following power-law is observed in the large scale region 

P{x) = Cx"^-^ for x>xo , (8) 

where xq is a certain threshold. This power-law is called Pareto's law [3] and the exponent /x 
is named Pareto index. Notice that Pareto's law does not hold below the threshold xq. The 
purpose of this section is to exhibit the distribution in the middle scale region under the detailed 
balance (7). 

In order to identify a statistical dependence in the growth rate of the past value, we examine 
the probability density faction of the profits growth rate in 2003-2004 firstly. In the database, we 
divide the range of xi into logarithmically equal bins as xi G 4 x [lO^^^-^^""^) , 10^+'^-^"] thousand 
yen with n = 1, 2, • • • , 20. Figure 5 shows the probability density functions for r = log^g R in 
the case of n = 1, • • • , 5, n = 6, • • • , 10, n = 11, • • • , 15 and n = 16, • • • , 20. The number of the 
firms in Fig. 5 is "22,005", "89,507", "85,020" and "24,203", respectively. 

From Fig. 5, we approximate logj^g ^(^l^i) by linear functions of r: 

logiQq{r\xi) = c(xi) — r for r>0, (9) 

logiQq{r\xi) = c{xi) + t-{xi) r for r<0. (10) 

Here q{r\xi) is the probability density function for r, which is related to that for R by Q{R\xi) = 
logiQq{r\xi) — r — logio(lnlO). These approximations (9)-(10) are expressed as tent-shaped 
exponential forms as follows 

Q{R\xi) = d{xi) for R>1 , (11) 

Q{R\xi) = d{xi) for R<1 , (12) 
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where d(xi) = lO'^'-^^^/ln 10. Furthermore, figure 5 shows that the dependenee of c{xi) on xi is 
negligible for n = 9, • • • , 20 . We assess the validity of these approximations against the results. 

Figure 6 represents the dependence of t±{xi) on the lower bound of eaeh bin .xi = 4 x 
]^Qi+o.2(n-i)_ YoT n = 17, • • • ,20, t±{xi) hardly responds to xi. This means that Gibrat's law 
(5) holds only in the large scale region of profits (x > xq)- ^ In contrast, t+(xi) linearly increases 
and t_(xi) linearly decreases symmetrically with log^gXi for n = 9, • • • , 13 . From Fig. 6, the 
slops are described as [13] 

t±{xi) = t±{xo) ± a In — . (13) 

Xq 

The parameters are estimated as follows: a ~ for xi > xq, a ~ 0.14 for Xjain < xi < xq, 
Xq = Ax 101+0-2(17-1) ^ 63,000 thousand yen and x^in = 4 x loi+0-2(9-i) 1,600 thousand 
yen. Notice that approximations (9)-(10) uniquely fix the expression of t-|-(.Ti) under the detailed 
balance [9]. This derivation is included in the proof in the next section. We call Eqs.(ll)-(13) 
Non-Gibrat's law. 

For n = 9, • • • , 20, the dependence of d{xi) on xi is negligible. In this case, the Non-Gibrat's 
law determines the probability density function of profits as follows 

P(xi) = Cxr^'^+i) e"'"'''' ^ for xi > x^in , (14) 

where t+(xo) — t_(xo) ~ /U [14]. This is the power-law in the large scale region (xi > xq) and 
the log-normal distribution in the middle scale region (xmin < xi < xq). The relations between 
parameters a, x in Eq. (1) and a, /x, xq are given by 

a = —^, /x = ^ln— . (15) 

Figure 7 shows that the distribution (14) fits with the empirical data consistently. Notice that 
the distribution cannot fit with the empirical data, if a is different from the value estimated in 
Fig. 6 (a = 0.10 or a = 0.20 for instance). 

These empirical data analyses are not restricted in the single term 2003-2004. Figures 8- 
10 show similar empirical data analyses in 2004-2005. The number of the firms in Fig. 8 is 
"20,669", "89,064", "88,651" and "26,887", respectively. Furthermore, empirical data analyses 
in 2003-2005 are shown in Figs. 11-13. The number of the firms in Fig. 11 is "18, 101", "76, 050", 
"75,451" and "22,456", respectively. In Figs. 9, 12, the parameters a, xq and Xmin in 2004-2005 
or 2003-2005 are estimated to be same values in 2003-2004. The distribution (14) in 2004-2005 
or 2003-2005 also fits with the empirical data consistently (Figs. 10, 13). 

^Fujiwara ct al. [10] prove that Parcto's law (8) is derived from the Gibrat's law (5) valid only in the large 
scale region and the detailed balance (7). In the derivation, linear approximations (9)-(10) need not be assumed. 
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3 Quasi-static log-normal distribution under the detailed quasi- 
balance 



In the previous section, the log-normal distribution in the middle scale region is exhibited by 
the Non-Gibrat's law under the detailed balance. The resultant profits distribution is empirically 
confirmed in data analyses of Japanese firms in 2003-2004, 2004-2005 and 2003-2005. The 
profits distribution (14) is static, because the derivation is based on the detailed balance (7) 
which is static time-reversal symmetry (Fig. 1-3). 

On the other hand, we have derived Pareto's law with annually varying Pareto index under 
the detailed quasi-balance [11]: 

/ \l/0 

Pi2ixi,X2) = Pi2{[^] ,ax/). (16) 

It is assumed that, in the quasi-static system, the joint probability density function has "a Xi^ <r^ 
X2 symmetry where ^ is a slope of a regression line: 

logio X2 = Q logio x\ -I- logio « ■ (17) 

The detailed balance (7) has the special symmetry Q = a = \. Because the detailed quasi-balance 
(16) is imposed on the system, Q is related to a as follows: 

= l-^logioa. (18) 

Here lO'" is a sufficient large value compared to the upper bound in which d and a are estimated. 

In Ref. [11], these results have been empirically confirmed by employing data on the assessed 
value of land [15]-[16] in 1983-2005. In the derivation, we have used the Gibrat's law (5) valid 
only in the large scale region without linear approximations (9)-(10). The purpose of this section 
is to show that the approximations uniquely fix a Non-Gibrat's law under the detailed quasi- 
balance. After that, we identify the quasi-static distribution not only in the large scale region 
but also in the middle scale region. 

By using the relation P\2{x\,X2)d.x\dx2 = Pir{xi, R)dxidR under the change of variables 
{xi,X2) {xi,R), these two joint probability density functions are related to each other 

Pir{xi,R) = a Pi2(xi,X2) , (19) 

where wc use a modified growth rate R = X2/{a xi^) . From this relation, the detailed quasi- 
balance (16) is rewritten as 

PiR{xi,R) = R-'PiR{(^^y\R-^) . (20) 

Substituting Pir{xi,R) for Q{R\xi) defined in Eq. (6), the detailed quasi-balance is reduced to 
be 

P(xi) _ 1 Q{R-^\{x2/af') 



P{{x2/af'') R Q{R\xi) 



(21) 
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In the quasi-static system, we also assume that Q{R\xi) follows the tent-shaped exponential 
forms (11)-(12). Under the approximations, the detailed quasi-balance is expressed as 



j^t+{xi)-t-{{x2/a)'^')+l 



(22) 



for R > 1. Here we use the notation P{x) = P{x)dix). By expanding Eq.(22) around = 1, 
the following differential equation is obtained: 

e \t+{x) - t-{x) + ll P{x) + xP '{x) = , (23) 

where x denotes xi. The same equation is obtained for < 1. The solution is expressed as 

P{x) = Cx-^ e-^(^) , (24) 

where 9 t+{x) — t-{x) = g{x) and / dx g{x)/x = G{x) . 

In order to make the solution (24) around R = 1 satisfy Eq. (22), the following equation 
must be valid for all R: 



-G{x) + G{R^'^x) = [t+{x) - t-{R^'^x) 
The derivative of this equation with respect to x is 

't+'{x)-R^" t.'{R^I' 



InR . 



g{x) ^ g{R'/^x) 



X 



InR . 



X X 

By expanding Eq. (26) around R=l, the following differential equations are obtained: 



X 



t+ {x) + t- {x) + t+ (x) + t_ (x) = , 



(1 - 361) t+ {x) + (2 - 361) t_ (x) + 3(1 - 9) x t+ (x) + 3(2 - 9) x t- (x) 



+ X 

The solutions are given by 
t+{x) = 
t.{x) = 



t+^^\x) + 2tJ'''>{x 



(3), 



. 



In^ X + (C+i - C_i) Inx + (C+o - C_o) , 
%^ln2x + C_ilnx + C_o . 



(25) 

(26) 

(27) 

(28) 

(29) 
(30) 



In order to make these solutions satisfy Eq. (25), the coefficients C_2 and C+i must be 
As a result, t±{x) under the detailed quasi-balance is uniquely fixed as 



X 



t±(x) = t±(xo) ± a In — . 

Xo 



(31) 



This is the same expression under the detailed balance. 

By using the Non-Gibrat's law (31), probability density functions Pi(xi), P2{x2) are uniquely 
reduced to be 



Pi(xi) = Ci xr''!-^ exp 
P2{x2) = C2 X2-'''-^ exp 



—9 a In 



2 ^ 

Xo 



-9 aln^ 



{x2/ay 

Xo 



(32) 
(33) 
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with 



= 9 



1^2 + 1 

in the large or middle scale region where the dependence of d{xi) on xi is negligible. 
Here we consider two log-normal distributions in the middle scale region: 



(34) 



-PlN2 



1 



1 



exp 



; exp 



hi^ (.ri/.xi) 

2(Tl2 

In^ (X2/X2) 



By comparing Eqs. (32)-(33) to (35)-(36), we identify 

1 



2(72^ 



9 a 

a. 



2(7i2 ' 

1 



— ^In — 

(71^ XX 

e 



1 axo 

Consequently, the relation between ai, 02 and is expressed as 



(35) 
(36) 

(37) 

(38) 

(39) 



This is the equation which quasi-static log-normal distributions satisfy. 



4 Conclusion 

In this study, by employing profits data of Japanese firms, we have exhibited the static 
log-normal distribution of profits in the middle scale region from a Non-Gibrat's law under 
the detailed balance. In the derivation, we have adopted two approximations observed in the 
database. One is that the probability density function of the profits growth rate is described 
as a tent-shaped exponential function. The other is that the value of the origin of the growth 
rate distribution divided into bins is constant. These approximations are confirmed in 2003- 
2004, 2004-2005 and 2003-2005. The resultant profits distribution fits with the empirical data 
consistently. This guarantees the validity of the approximations. 

This static discussion has been applied to a quasi-static system. We have derived a quasi- 
static log-normal distribution in the middle scale region from a Non-Gibrat's law under the 
detailed quasi-balance. In the derivation, we have also assumed two approximations confirmed 
in the database. Even in the quasi-static system, the approximations uniquely fix the Non- 
Gibrat's law to be the same expression under the detailed balance. The resultant distribution 
is power-law with varying Pareto index in the large scale region and the quasi-static log-normal 
distribution in the middle scale region. Notice that not only the change of Pareto index \x but 
also the change of the variance of the log-normal distribution a depends on a parameter of the 
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detailed quasi-balance 9. In other words, /x and a are related to each other. This phenomenon 
is suggested in Ref. [17]. 

First issue in the future is to confirm this quasi-static derivation by data analyses. For this 
aim, long-term economic data must be investigated in the middle scale region where dynamical 
transitions arc observed. We believe that the analytic results in this study are confirmed in such 
a database, if two approximations assumed in this paper are applicable. 
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Figure 1: The scatter plot of all firms in the Figure 2: The scatter plot of all firms in the 

database, the profits of which in 2003 {xi) and database, the profits of which in 2004 (xi) and 

2004 (^2) exceeded 0, xi > and X2 > 0. The 2005 (^2) exceeded 0, xi > and X2 > 0. The 

number of the firms is "227,132". number of the firms is "232,497 
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Figure 3: The scatter plot of all firms in the Figure 4: Probability density functions of prof- 
database, the profits of which in 2003 (xi) and its in 2003-2005. In the large scale region, 
2005 (X2) exceeded 0, xi > and X2 > 0. The Pareto's law is observed. Pareto index is es- 
number of the firms is "197,867 timated to be nearly 1. 



10 



-2-1012 

TheLogarithmofProfitsGrowthRater=log ,oX A , 



-2-1012 

TheLogarithmofProfitsGrowthRater=log ,oX A , 





Figure 6: The relation between the lower bound 
of each bin xi and t±{xi) with respect to 
the profits growth rate from 2003 to 2004. 
Prom the left, each data point represents n = 
l,2,---,20. 



Figure 7: The probability distribution function 
of profits P{xi) for firms, the profits of which 
in 2003 (xi) and 2004 (X2) exceeded 0, xi > 
and X2 > 0. 
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Figure 9: The relation between the lower bound 
of each bin xi and t±{xi) with respect to 
the profits growth rate from 2004 to 2005. 
Prom the left, each data point represents n = 
l,2,---,20. 



Figure 10: The probability distribution func- 
tion of profits P{xi) for firms, tlie profits of 
which in 2004 (xi) and 2005 {x2) exceeded 0, 
> and X2 > 0. 
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Figure 12: The relation between the lower 
bound of each bin xi and t±{xi) with respect 
to the profits growth rate from 2003 to 2005. 
From the left, each data point represents n = 
l,2,---,20. 



Figure 13: The probability distribution func- 
tion of profits P{xi) for firms, tlie profits of 
which in 2003 (xi) and 2005 (X2) exceeded 0, 
xi > and X2 > 0. 
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